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Is the closest product state to a symmetric entangled multiparticle state also symmetric? This 
question has appeared in the recent literature concerning the geometric measure of entanglement. 
First, we show that a positive answer can be derived from results concerning symmetric multilinear 
forms and homogeneous polynomials, implying that the closest product state can be chosen to be 
symmetric. We then prove the stronger result that the closest product state to any symmetric 
multiparticle quantum state is necessarily symmetric. Moreover, we discuss generalizations of our 
result and the case of translationally invariant states, which can occur in spin models. 
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I. INTRODUCTION 

Entanglement is a key phenomenon in quantum me- 
chanics and its quantification is vital for the field of quan- 
tum information theory. Many entanglement measures 
have been proposed for the two-particle as well as for 
the multiparticle case [H, g. Virtually all of the pro- 
posed entanglement measures, however, suffer a serious 
drawback: They are very difficult to compute as their 
definition contains optimizations over certain states or 
quantum information protocols 3]. Such optimizations 
can be performed successfully for special cases only, for 
instance, if the density matrix under investigation pos- 
sesses a high symmetry or belongs to a special family, 
e.g., with low rank [1]. 

An often-used entanglement measure for multiparticle 
systems is the geometric measure of entanglement 
For a given multiparticle state li/j), one first considers the 
closest fully separable state {(p) = \a)\b)\c) ■ ■ ■ in terms of 
the overlap 



max 

) = \a)\b)\c)- 



(1) 



and then defines the geometric measure of the pure state 
as 



(2) 



Sometimes, the geometric measure for pure states is also 
taken as ecH^j)) = — 21og2 G'(?A). Based on this defini- 
tion, the geometric measure is extended to mixed states 
via the convex roof construction: For a given density ma- 
trix g one minimizes over all possible decompositions of 
g into pure states g = J^kPi' \ 'i'k){(t>k\, where the pk form 
a probability distribution. 



Ecig) = min S2pkEG{\(l>k))- 

Pk,\4>k) , 

k 



(3) 



Clearly, also this optimization is not straightforward to 
compute. 



The geometric measure has become one of the widely 
used entanglement measures for the multiparticle case. 
It fulfills all the desired properties of an entanglement 
monotone 0]. Moreover, it has a physical interpretation 
of quantifying the difficulty in distinguishing multiparti- 
cle quantum states by local means [6,]. It has also been 
used to study quantum phase transitions in spin mod- 
els 0, HI and the usefulness of states as resources for mea- 
surement based quantum computation Q. The value of 
Eg has been computed for many pure states [13, [HI, [I3l , 
and the convex roof for some important cases has been 
calculated in Refs. [slfTst. 

If one considers the optimization problem in Eq. ([T]), 
a natural question arises whether for a symmetric state 
\tp) the closest product state can be chosen symmetric, 
i.e., 10) = |a)|a)|a) • • • . If this is true, it drastically sim- 
plifies the calculation of the geometric measure for pure 
symmetric states, as the number of parameters in this 
optimization then does not depend on the number of par- 
ticles anymore. Recently this problem drew considerable 
attention in quantum information theory and some ef- 
fort was made to prove it. For example, it has been 
used as a conjecture in Ref. [1]. In Ref. [IH it has been 
proved for two particles that there is always a symmetric 
state which gives the maximum value (but it can hap- 
pen that also non-symmetric states yield the same value) 
and a first attempt for the A^-particle case was given. 
Quite recently, special cases of this conjecture have been 
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14l |. and related conjectures have been formu- 



verified 
lated 

In this paper, we investigate the conjecture from 
several perspectives. We show that a result on A^- 
homogeneous polynomials over Banach spaces can be ap- 
plied to the above problem and proves that the maximum 
in Eq. ([T]) can be achieved by a symmetric state. How- 
ever, this result does not allow to conclude that only sym- 
metric solutions exist. We then go on to show that the 
optimal state maximizing G^tp) is necessarily symmetric 
for three or more particles. Finally, we will discuss con- 
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sequences and generalizations of our results, concerning, 
among others, the maximization of the expectation value 
of symmetric positive operators. 

II. THE MAIN RESULT 

In this section we will first apply a result from the the- 
ory of homogeneous polynomials to the conjecture, prov- 
ing that a symmetric state attains the maximum. Then 
we proof a stronger version, stating that the maximizing 
state is necessarily symmetric. Let us fix our notation. 
We write the overlap with a product state as an evalua- 
tion of a corresponding A'^-linear form, i.e., 

{ip\ai,a2, . . . , aAf) =: ■>jj{ai,a2, . . . , aw). (4) 

A symmetric state {ip) then corresponds to a symmetric 
TV- linear form ip. For a fc- level system, the vectors 
are elements of (D*^ and we will always assume they are 
normalized. Furthermore, the expression a = P denotes 
equality up to a phase, i.e., there is a phase ip such that 
a = e*^/3. 

For any symmetric A'^-linear form there exists an as- 
sociated iV-homogeneous polynomial tp [16] by virtue of 
the mapping 

ip: a 1-^ ip{a, . . . , a). (5) 
The norm of this polynomial is defined by 

W^PW = max\iip{a)\. (6) 

a 

In this context the quantity G{ip) we introduced above 
is denoted by ||^||. 

The relation between the norms HV^H and HV'II is stud- 
ied in the theory of polynomials over Banach spaces, cf. 
Ref. [13] ■ For a Banach space £ the polarization constant 
c{N,£) is the smallest positive number such that 

U\\<c{N,£)m (7) 

holds for every symmetric iV-linear form ip over £. 

For finite dimensional real and complex Hilbert spaces 
this polarization constant is known to be 1 , cf. Ref. llSl . 
ITgj . An exhaustive discussion can be found in Ref. [l7j |. 
In our particular case £ = (D*^ it follows that symmetric 
product states maximize the overlap with any symmetric 
state, but leaves open if also non-symmetric states attain 
the maximal value. This is in fact never the case for 
A'' > 3, as we summarize in 

Lemma 1. Let rp ^ be a symmetric N -linear form over 
(D*^ with N > 3 and let the vectors ai, . . . ,aN maximize 
\ip\, i.e., 

G(V') = |'(/'(Q;i,a2, . . . ,aAr)|. (8) 

Then the vectors ak are equal up to a phase, in other 
words, the span of ai, ajsr is one-dimensional. 



In order to prove Lemma [U we will first consider the 
situation where N = 2. The following Lemma and its 
proof were already given Ref. [Tl| . however, the proof 
provides some observations that are essential in order to 
establish our main result. 

Lemma 2 For any symmetric two-linear form tp 

over (D*^ we can find a vector a such that 

G(V) = |^(a,a)|. (9) 

Rephrasing the statement of the Lemma, when maxi- 
mizing \^\ for two particles, the maximum can by reached 
by a symmetric choice of vectors, although a solution 
maximizing ji/'l is not necessarily symmetric. 

Proof of Lemma\^ In a fixed orthonormal basis {6^} the 
symmetric quadratic form -0 is represented by a symmet- 
ric matrix '^ij := ip{hi,hj). Then we have '>p{a,pi) = 
a^'^ p}, where on the right-hand side a and /3 are column 
vectors with coefficients in the basis {hi}. For a complex 
symmetric matrix ^ = Takagi's factorization theo- 
rem states that can be written as 

* = U'^DU, (10) 

with a unitary matrix U and a diagonal matrix D = 
diag(ri, . . . ,rfc) where the non- negative values are in 
decreasing order, ri > r2 > • • • > rj, > 0. In this form, 
a symmetric choice of a and /3 maximizing \ip\ becomes 
evident, namely a = pi = Wei where ei = (1, 0, . . . , 0)"^. 
Hence G'('0) = r-i. □ 

Let us make some remarks on the remaining freedom 
in the choice of a and /3. First, we note that if ri > r2 
the only choice to reach the maximum ri is the one given 
in the proof above. Hence for this case the maximizing 
solution is unique (up to a phase) and symmetric. 

Otherwise, consider the case ri — r2 = ■ ■ ■ = 
rd. We then say ip is degenerate and define Ri 
spanj.({ei, . . . ,6^}). If G{'tp) = \^p{a,p))\, we can always 
write a = U^e* and p) = We with some e G Ri, where 
the vector e* e Ri denotes the vector obtained from e 
by complex conjugation in the given basis. The case that 
e = e* then corresponds to the symmetric solutions. 

The symmetric maximizing solutions therefore corre- 
spond via U to real vectors in i?i (up to a phase), and 
having a non-symmetric solution a ^ /3 of the maximiza- 
tion implies degeneracy. Moreover, in case of degeneracy 
we find a continuum of inequivalent asymmetric as well 
as symmetric solutions. The following observation will 
be needed in the main proof and expresses this fact. 

Observation 3. Given a symmetric two-linear form tp 
with 

G(V') = |V(a,/3)|, (11) 

where pi ^ a, we can always find two orthonormal vectors 
5i and 62, such that: 
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(i) Si and 62 span the same space as a and (3, in par- 
ticular a, P G span^({Si, 62})- 

(11) We have G{i(;) = \i^{di,5i)\ = \ip{S2,S2)\ = 
\-ip{T],T])\ = \ilj{^i,fi')\, where 

V (51 + 52)/V2, 

ti:^{5i + i62)/V2, (12) 
// := {Si ~ iS2)/V2- 

(Hi) The vectors Si,S2, rj, /i' do not equal a, even modulo 
a phase. 

Proof. Consider G(?/') = |V'(a,/3)| = \^{We*,We)\ with 
the unitary matrix U from Takagi's factorization of 5*, 
where e ^ e* as ^ a by assumption. We can 
choose two real orthonormal vectors /i , /2 such that 
span((;,({/i, /2}) = span<p({e, e*}). They can be obtained 
from the real and imaginary parts of e and e* [2l| . Hence 

5i = U^fi and S2 = [/t/2 (13) 

is a vahd choice of orthonormal vectors fulfiUing (i), each 
providing a symmetric maximization of \Tp\. The vectors 
Si,S2 as weU as the vectors r/,fi,^' derived from Si,S2 
fulfill (ii) by construction, according to the observations 
above and because 77 = U^fi + /2)/\/2, where (/i + 
/2)/a/2 € Ri (likewise for /i, /i'). Finally, we have enough 
freedom to choose /i, /2 to satisfy (iii). □ 

After these preliminaries, we are ready to prove 
Lemma [TJ 

Proof of LemmaUl The proof consists of two parts. In 
Part 1 we prove the case iV = 3. In Part II we extend the 
result to arbitrary > 3. 

Part I. Assume a maximizing set of vectors {a,/3,"f} 
has been found, 

G(^) = |V(a,/3,7)l- (14) 

We show that the assumption dim[span£,({a, 7})] 7^ 1 
leads to V' = 0. 

Without losing generality, we hence assume that 7^/3. 
Then we have a degenerate quadratic form ip^a, •, ••). Us- 
ing Lemma [5] we obtain a symmetric maximizing solu- 
tion a ^ a (due to 7 ^ /?, u is not unique). Hence 
we have G{tp) = \^{a,a,a)\. From Observation [3l ap- 
plied to the quadratic form ?/'(•,••, cr), we take the vectors 
Si,S2,ri, n, fi' with the properties as stated. With these 
vectors we define the 2 x 2-matrices A,B,N,M via 

Aki := 'ip{Si,Sk,Si), 

Bki ■■= ip{S2,Sk,Si), 

Nki ^(77, 4, Si) = {Aki + Bki)/V2, ^^^^ 
Mm ?A(m, 4, Si) = {Aki + iBki)/V2. 

As (7, (5i, (52, ?7, and fi' are in spau;]-. ({(5i, 4}) the ma- 
trices A, B, N, M correspond to two-forms which as- 
sume the maximum ri = G{^p) on this span. Since 



Si,S2,ri, n' ^ (T, the quadratic forms A,B,N,M are 
degenerate with the value ri. In terms of the Tak- 
agi factorization, we can write A — U'^DUj^ and B = 
UgDUg, etc. with D = diag(ri,ri). It follows that 
A^A = B^B ^ N^N = M^M = D^, which implies 
A^B + B'^A = i{A^B - St A) = 0. Hence B^ A = 0, 
which is only possible for ^ = 0. 

Part II. The extension to the case A^ > 3 is proved 
as follows. If G{'4!) = |^/'(ai, q;2, . . . , ckAf)! we define a 
symmetric 3-linear form by 

V'^ := -0(ai,a2, • ■ • ,afc, ■ • ■ ,aAr) (16) 

where denotes omission, i.e., ai,Q;2,afc are omitted. 
As in Part 1, we have G{'ip) = |V'^.(q!i, a2, afc)| and thus 
ai = a2 = Oik for all k. Hence all vectors ai, . . . ,aN 
must be the same up to a phase. □ 

Lemma [T] is stated for A^-linear forms over a complex 
vector space. If is real, one can also consider the 
maximization over real product vectors. In general this 
will yield a different result than the complex case [2^ . 
Since the polarization constant is 1 also in the real case, 
one can find a a symmetric state among the real prod- 
uct states which attains the maximum. In contrast to 
the complex case, the maximizing state is, however, not 
necessarily symmetric for three particles. A counterex- 
ample is IV") = (|001)-|-|010)-f|100)-|lll))/2 where the 
maximum of 1/2 is also attained by — |001). 

The following Lemma provides an additional very sim- 
ple proof that c{N, H) — 1 for A^ — 2^, and is based on a 
symmetrization procedure. 

Lemma 4. Let ip he a symmetric N-linear form with 
N ^ 2^ and let the vectors ai,...,Q!jv maximize {tpl. 
Then there exists a normalized vector C in the span of 
Q!i, . . . , QfAT such that 

|V(«i,...,«jv)| = |^(C,...,C)I- (17) 

This statement holds for real and complex Hilbert spaces. 

Proof. Let (•,••) denote the scalar product and we 
define the constant A = ^p{ai, a2, a^, . . .). Then 
Tp{-,a2,a3, . . .) = A(ai,-) and ?/'(ai, ' , "3, ■ • ■ ) = 
A(q!2, •) (cf. Eq. (6) in Ref. Q). Using the symmetry 
and linearity of ip we arrive at 

V'(-,/3i,a3,...) = A(/3i, •) (18) 

and hence 

A = V(/3i,/?i,a3,...), (19) 

where /3i = {ai + a2)/\\ai -I- Q;2||- (If ai = —c(2, we set 
/3i = ai and replace A by —A.) 

We now repeat this procedure, first yielding A = 
V'(/3i,/?i,/33,/33, • ■ • ) and then A = -0(71,71,71,71, . . . ), 
where P3 is defined analogously to f3i and 71 — {f3i + 
P3)/\\f3i + PaW- In the second step we applied the sym- 
metrization to the first and third argument as well as 
to the second and fourth argument. Since N = 2^, 
we can complete this symmetrization and arrive at A = 
ipiC, • ■ ■ ) C) foi' some C, in the span of ai, . . . ,aN. □ 
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III. DISCUSSION 



A. Physical interpretation of the proof 



be handled with care. For some models, it seems to be 
true Isl. [23l . for other models (like the Majumdar-Ghosh 
model [2^) one can directly check that it is wrong. 



An interpretation of the proof of Lemma [T] in physical 
terms is the following. The matrices A, B, N and M in 
Eq. are representations of the state after one site 
has been measured out and the remaining state has been 
projected onto a two-dimensional subspace. The values 
ri correspond to Schmidt coefficients of this remaining 
state and, as they are equal, the state corresponds to a 
Bell state. The proof of Lemma [1] shows that for qubits 
it is impossible to create a state of three particles that 
is both symmetric and always results in a Bell-pair like 
state after an arbitrary measurement on one site. 



B. Translationally invariant states 

It is interesting to ask whether also for translationally 
invariant states the maximum is attained in a symmetric 
state, as such states occur naturally in the analysis of 
spin models. This has sometimes been assumed when 
investigating the geometric measure in condensed matter 
systems. 

First, a counterexample for this conjecture is the state 
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1^-) = ^(10101) + 11010)) 



(20) 



for which the closest separable states are the non- 
symmetric states 1 0101) and |1010). In fact, one can 
find translationally invariant states which are orthogo- 
nal to any symmetric product state, e.g., ^ (|0101) — 
|0011) + all translations). 

This situation gets worse as the number of particles 
increases. Let T denote the subspace of translationally 
invariant states for N qukits and let 5 C T be the per- 
mutationally symmetric subspace. Then any state in 
X — T n S-^ - the orthocomplement of 5 in T - has 
a vanishing overlap with any symmetric product state, 
hence the closest product state is not symmetric. The 
dimension of T is given by [23] 



dim(r) = l^^(j)fc^/^ 



(21) 



where (p denotes Euler's totient function and the sum- 
mation is over all divisors j of N. For S we have 



dim(5) ^ [ 



(22) 



Therefore, if fc then the dimension of the subspace 
X is roughly given by {k^ — N'^)/N and the fraction of 
states where the conjecture holds shrinks rapidly as the 
number of particles increases. 

Concerning the analysis of entanglement in spin mod- 
els, this shows that the assumption that the closest sep- 
arable state to the ground state is symmetric, has to 



C. Operators of higher rank 

We now consider generalizations of our results. Let 
Tig be the projector onto the symmetric subspace S. An 
operator A is permutationally symmetric if it acts on the 
symmetric subspace only, i.e., it fulfills A — YlgAIls. 
A is called permutationally invariant if it is invariant 
under permutation of the particles (the latter is a weaker 
condition than the former [1^). We hence define for an 
observable X 



G(X):= max MX\ip)\ 

\ip) = \a)\b)\c)--- 

Gs{X):= naax K^^lAl^.)! 

\ip} = \a}---\a) 



(23) 



Such optimizations occur naturally in the construction of 
entanglement witnesses or in the estimation of entangle- 
ment measures via Legendre transforms 

To study the relation of these quantities, we can write 
{ip\X\ip) = Tr [X \(p){(p\] as an evaluation of a correspond- 
ing A^-linear form ^ over Sj'' (the Banach space of Her- 
mitian matrices of dimension k equipped with the trace 
norm) due to 



,AN)^Tr[XAi 



'A 



(24) 



Then any permutationally invariant operator X cor- 
responds to a symmetric A^-linear form ^ and an A^- 
homogeneous polynomial ^. We now define ||^|| and 
ll^ll analogously to Section |IT] (using \\A\\ = 1 as nor- 
malization condition). It is straightforward to see that 
G{X) = ll^ll and GsiX) =Jji||- As the polarization 
constant can be shown to be [271 



c{N,Sj'') = N^/m for N <k, 



(25) 



the quotient G{X) /Gs{X) can get arbitrarily large as N 
and k increase. 

At the end of this section, we will provide a further 
explicit example. Let us first discuss some cases where 
symmetry assumptions do hold: 

Corollary 5. (i) If X is a positive permutationally sym- 
metric observable then G{X) can be attained by a sym- 
metric state. 

(a) If X is a permutationally invariant N-qubit observ- 
able that contains only full correlation terms, then G{X) 
can be attained by a symmetric state. 

Proof, (i) We note that 

G(A)< max max (26) 

\i>) = \bi)---\br,) \ip) = \ai)---\a„) 

Fixing the (unnormalized) state A|-0) = n^AII^lT/;) 
is symmetric and by virtue of Lemma [T] the maximum 
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is reached by a symmetric state \a,...,a). Repeat- 
ing the reasoning with the fixed state \ip) = |a, . . . , a), 
we get G{X) < max|b^ maxi^^ |(a, . . . , a|X|6, . . . , The 
fact that for positive operators 2|(q!|P|/3)| < (q;|P|q;) + 
{f3\P\(3) < 2max{(a|P|a), (/3|P|/3)} holds for arbitrary 
|a) and |/3) proves that G{X) < GsiX), hence G{X) = 
Gs{X). 

(ii) An iV-qubit operator X contains only full correla- 
tion terms if X = J2i,j,...£{x,y,z} ® crj <Si ■ ■ ■ . Note 
that here ao — 1 does not occur; a physically relevant and 
well known example for such an operator X is the Mermin 
inequality. As X is permutationally invariant, A is equiv- 
alent to a symmetric A'^-linear form over R'^. The Bloch 
representation for qubits implies that here the maximiza- 
tion is equivalent to finding max^.gj^a |A(ri, r2, . . . , rAr)| 
where the vectors are the corresponding Bloch vectors. 
Since the polarization constant for real Hilbert spaces is 
1 [13, m, 111, the assertion follows. □ 

Let us conclude with some examples where symmetry 
assumptions do not hold. If X is symmetric but not 
positive, then the maximum G{X) is, in general, not 
attained by a symmetric state. A counterexample for 
two qubits is a: = 6 - |00)(00| - 2 |11)(11| with 

1^^+) = (|01) + \10))/V2. Then G{X) = 3 (we can take 
|0o) ~ |01)) while the maximum for symmetric prod- 
uct states is 34/15. Also, if X is permutationally in- 
variant (and even positive) then the maximum G{X) 
is, in general, not attained by a symmetric state. A 
counterexample is the singlet state, X = \ip~){ip~\ with 
= (|01) — \10))/^/2. This operator is invariant under 
permutation of the particles, but it docs not act on the 



symmetric space. It has G{X) = 1/2, but restriction to 
the symmetric would yield again G{X) = 0. This 
clarifies some questions raised in Ref. ^5,] . 



IV. CONCLUSION 

In conclusion, we have discussed a widely used con- 
jecture concerning the geometric measure of entangle- 
ment. Our results not only simplify the calculation of 
the geometric measure for symmetric states, but they 
also have applications to subjects in condensed matter 
physics. Furthermore, from a mathematical perspective, 
the quantity G, as defined in the introduction, is known 
as the injective tensor norm |28|. On the one hand, this 
norm is of central importance in tensor analysis, since 
there are, as Grothendieck showed, fourteen inequivalent 
natural tensor norms derivable from G [l^ . On the other 
hand, this norm has also occurred in the discussion of the 
maximal output purity of quantum channels lid\. So we 
believe that the study of tensor norms can yield further 
interesting insights in quantum information theory. 
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